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Abstract —  We  present  some  2D  and  3D  texture  fea¬ 
tures  computed  from  the  grey  values  of  MRI-T1  data. 
The  features  show  strong  correlation  with  the  score  in 
the  Mini  Mental  State  Examination  (MMSE)  used  rou¬ 
tinely  to  help  diagnose  Alzheimer’s  disease. 

I.  Introduction 

In  this  paper  we  use  the  method  of  graylevel  dependence 
histograms  (GLDH)  as  defined  by  Chetverikov  for  2D  [2] 
and  generalized  to  3D  by  Kovalev  and  Petrou  [9]  and 
derive  texture  anisotropy  features  from  MRI  data,  that 
correlate  well  with  the  result  of  Mini  Mental  State  Exam¬ 
ination  (MMSE)  which  is  routinely  used  to  help  diagnose 
Alzheimer's  disease  (AD)  [3]. 

Freeborough  and  Fox  [5]  used  the  2D-GLDFI  method 
applied  to  MRI  data,  but  the  features  they  extracted  did 
not  measure  anisotropy.  Their  MRI  scans  were  acquired 
from  40  normal  controls  and  24  AD  patients  and  a  clas¬ 
sification  rate  of  91  %  was  obtained.  Mathias  et  al.  [11] 
also  used  the  2D-GLDH  method  applied  to  spinal  cord 
MR  images.  They  found  significant  differences  in  texture 
between  normal  controls  and  multiple  sclerosis  patients, 
and  also  a  significant  correlation  between  texture  and  dis¬ 
ability. 

For  this  study  we  used  25  scans  MRI-T1  weighted,  with 
coronal  orientation.  Each  data  cube  is  180  x  180  x  124 
pixels  and  the  voxel  size  in  mm  is  (0.9375,0.9375,1.5). 
Thirteen  of  these  scans  are  Alzheimer’s  patients  identified 
as  ADI,  AD2,  ...,  AD13  and  the  other  twelve  are  control 
volunteers  identified  as  C01,  C02,...,  C012.  The  control 
group  was  matched  in  age  and  gender  with  the  AD  group. 
The  mean  age  of  the  patients  with  AD  at  the  time  of 
their  scanning  was  56.77  with  an  age  range  from  39  to 
72.  While  the  mean  age  of  the  control  volunteers  was 
58.33  with  an  age  range  of  47  to  72. 

The  Mini  Mental  State  Examination  (MMSE)  score 
for  each  of  the  patients  is  also  available.  The  MMSE 
is  a  test  used  to  detect  dementia,  although  it  is  not 
specific  to  AD.  The  maximum  score  is  30  and  controls 
will  typically  score  29  or  30/30.  Scores  between  10  and 
24  are  considered  mild  to  moderate  dementia  cases,  and 
scores  below  10  show  severe  dementia.  Table  I  displays 
the  scores  for  each  patient.  Two  of  the  scores  do  not 
match  with  the  clinical  diagnosis:  AD3  and  C02.  The 
AD3  scan  is  interesting  as  it  comes  from  a  patient  who 
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TABLE  I 

MMSE  Scores  for  each  patient 


was  imaged  just  before  the  onset  of  the  first  clinical 
symptoms,  and  at  a  time  when  there  might  have  been 
ongoing  structural  brain  changes  -  hence  the  high  score 
for  that  subject. 

The  scans  were  segmented  to  isolate  the  brain  from 
external  structures  (eyes,  ventricles,  bones,  etc.)  [6]. 
The  brains  were  further  segmented  to  isolate  the  white 
and  grey  matter  as  well  as  the  border  between  these 
two  types  of  tissue  [1] . 

The  texture  analysis  technique  we  use  effectively 
counts  the  number  of  pairs  of  voxels  that  appear  in 
the  same  relative  position  and  have  certain  fixed  grey 
values.  Therefore  the  relative  grey  values  of  the  vox¬ 
els  are  extremely  important.  In  that  respect  MRI  brain 
images  may  not  be  compatible  with  each  other,  depend¬ 
ing  on  the  initial  conditions  of  the  scanning  process.  A 
process  of  normalization  is  used  in  order  to  have  the 
same  relative  grey-level  values  for  different  scans.  The 
smallest  grey-level  value  for  the  segmented  scan  is  as¬ 
signed  1  and  the  highest  255,  and  0  labels  the  voxels 
which  do  not  belong  to  the  region  of  interest. 

II.  3D  Texture  Representation 

First  we  start  by  defining  a  coordinate  system  in  our 
data  cube:  ( x,y )  axes  are  on  the  plane  of  each  scan 
slice,  ^  is  along  the  axis  perpendicular  to  each  slice. 
Azimuthal  angle  </>  is  measured  on  the  ( x ,  y )  plane  from 
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the  direction  of  the  x  axis.  A  pair  of  values  (<f>,  z )  defines 
a  unique  orientation  in  the  3 D  space. 

To  characterize  the  texture  in  the  data  we  compute 
the  quantity  [7]: 


255  255 

h(4 <>,  #  <0  =  E  E^fc  ~  l)2C(k,  l;  <j>,  z ,  d)  (1) 

k  1  1=1 

where  C(k,l;<f>,z,d )  is  the  number  of  pairs  of  voxels 
which  are  at  distance  d  from  each  other,  along  the  di¬ 
rection  defined  by  ((f),  z),  with  one  member  of  the  pair 
having  a  grey  value  k  and  the  other  l. 

If  the  data  are  completely  isotropic,  function 
h(<f>,  z;  d)  must  be  independent  of  direction  (<f>,  z)  and 
when  it  is  therefore  plotted  as  the  polar  radius  along 
direction  (<p\  z)  in  3 D  space,  it  must  represent  the  sur¬ 
face  of  a  sphere.  Any  deviation  from  the  shape  of  a 
sphere  indicates  anisotropy  in  the  data.  For  practi¬ 
cal  purposes,  function  h(<f),z;d )  is  sampled  at  a  dis¬ 
crete  number  of  points  (<f>i,  Zj)  where  i  =  1,2 , AT  and 
j  =  1,2, ...,  M.  The  sample  points  are  chosen  in  such 
a  way  that  each  represents  a  patch  of  constant  area 
on  the  surface  of  a  sphere  [10].  We  represent  the  dis¬ 
cretized  version  of  function  h(cf),  z:  d)  by  H(i,j )  and  we 
understand  that  it  is  calculated  for  various  values  of 
distance  d.  When  plotted  in  3 D  space  (for  fixed  d)  it 
represents  a  closed  digital  surface  which  we  call  indi- 
catrix.  An  example  of  such  an  indicatrix  is  shown  in 
figure  1. 


Alzheimer  Scan.  GLDH,  contrast,  d  =  0.9375  mm 
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Fig.  1.  Example  of  a  3D  indicatrix 

Our  purpose  is  to  measure  how  much  the  shape  of 
this  indicatrix  deviates  from  the  shape  of  a  sphere.  This 
will  be  a  measure  of  anisotropy  in  the  data.  In  the  next 
section  we  define  some  features  which  are  designed  to 
measure  exactly  that. 

III.  Feature  Extraction 

Three  features  proposed  by  Kovalev  and  Petrou  [9] 
are  used  to  analyse  the  shape  of  the  3 D  indicatrix. 

FI  is  the  Anisotropy  Coefficient 

Fl  =  (2) 

Fkmin 


F2  is  the  Integral  Anisotropy  Measure  or  standard 
deviation: 


Fo  = 


'IXi  Ef=1  ( H(i,j)-Hmy 
NM 


(3) 


F3  is  the  Local  Mean  Curvature  or  the  average  value 
of  the  Laplacian: 


f3  = 


UH(i,j+  1)) 


N(M  -  2) 


(4) 


where  Hmin  and  Hmax  are  the  minimum  and  max¬ 
imum  value  of  the  indicatrix  respectively,  and  Hm  is 
the  mean  value  of  all  points  of  the  indicatrix. 

Another  set  of  features  can  be  extracted  by  expand¬ 
ing  the  indicatrix  in  terms  of  spherical  harmonics  [8]. 
The  coefficients  of  such  an  expansion  can  characterise 
any  3D  closed  surface. 

The  equation  that  describes  the  spherical  harmonic 
expansion  of  a  surface  R(0,  cf>)  is  as  follows: 


L  l 

R(6A )  =  E  E  Ai,mUlim(0,<l>)  +  Bl,mvl,m(9,<f>)  (5) 

1=0  m= 0 

Ui,m(9,  (p)  =  cos(m^)P™(cos6»),  (6) 

Vi,m(0,  <f>)  =  sin(m</>)P/" (cos  9)  (7) 

where  9  is  the  latitude  angle,  (f>  is  the  longitude  angle, 
and  P™(cos<9)  is  the  associated  Legendre  polynomial 
[12]. 

Coefficient  Ao,o  is  the  mean  radius  of  the  indicatrix, 
i.e.  Ao,o  =  Hm.  All  the  other  non-zero  coefficients 
represent  different  types  of  anisotropy.  Figure  2  shows 
the  type  of  anisotropy  represented  by  some  low  order 
coefficients. 


IV.  Features  correlating  with  the  MMSE 
score 

We  extracted  anisotropy  features  from  four  differ¬ 
ent  brain  regions:  the  whole  brain,  white  matter,  grey 
matter,  and  the  border  between  white  and  grey  matter. 
We  define  the  border  between  grey  and  white  matter 
by  performing  a  dilation  of  3  x  3  x  3  in  the  voxels  from 
the  transition  white-grey  matter.  In  every  single  re¬ 
gion  five  different  distances  d  were  used:  0.9375mm1, 
1.5mm,  2 mm,  2.5mm  and  3mm,  ie  5  separate  indi- 
catrices  were  contracted  from  each  region.  Table  II 
presents  the  summary  of  the  experiments  performed 
and  the  best  features  identified.  The  last  column  gives 
the  correlation  coefficient  of  the  particular  feature  with 
the  MMSE  score. 

1Note  that  due  to  the  intraslice  resolution  being  0.9375  mm, 
the  use  of  d  =  0.9375  mm  corresponds  to  2D  analysis  if  not  grey 
level  interpolation  takes  place.  However,  since  we  perform  3D 
trilinear  interpolation  to  find  grey  values  at  arbitrary  positions, 
we  extract  3D  information  even  when  d  is  less  than  the  interslice 
distance. 


Spherical  Harmonic  A& 5  Spherical  Harmonic  B4  4  TABLE  II 

The  best  features  identified  by  all  experiments 
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Fig.  2.  Spherical  Harmonics  superimposed  on  the  unit  sphere 

The  features  which  best  correlate  with  the  MMSE 
score  were  computed  from  the  grey  matter  (shown  high- 
lated  in  table  II).  Figure  3  shows  feature  A±yi  plotted 
versus  the  MMSE  score,  where  the  correlation  between 
these  two  seemingly  independent  variables  is  very  clear. 


Component 

Best 

Features 

Correlation 

with  MMSE 

Whole  brain 

!  ^5,4 1 

d  = 

0.9375 

-0.749 

|  Z?3,2 1 

d  = 

2.5 

-0.747 

1^6,1 1 

d  = 

2.5 

-0.739 

White  matter 

1^45,3 1 

d  = 

2 

-0.750 

1-4-6,  si 

d  = 

0.9375 

-0.647 

Grey  matter 

iH 

I— 1 

77 

d  = 

0.9375 

-0.876 

l-As.sl 

d  = 

0.9375 

-0.845 

F3 

d  = 

0.9375 

-0.830 

f3 

d  = 

1.5 

-0.822 

iH 

iH 

7L 

d  = 

1.5 

-0.818 

iH 

rH 

27 

d  = 

3 

-0.818 

\A&M 

d  = 

0.9375 

-0.817 

|  -4-4,0 1 

d  = 

2 

-0.813 

f3 

d  = 

2 

-0.810 

Grey/white 

•Bi,i 

d  = 

2.5 

-0.734 

border 

Bi,i 

d  = 

3 

-0.729 

B1A 

d  = 

2 

-0.701 

lar,  we  are  examing.  1.  Whether  2D  are  adequate  or  3D 
data  are  necessary  for  the  derivation  of  such  features. 

2.  Whether  the  correlations  we  found  are  stronger  or 
weaker  than  correlations  with  the  age  of  the  subjects. 

3.  Whether  the  features  we  computed  reflect  more  the 
shape  of  the  region  of  interest  than  the  texture  inside 
it.  Results  of  the  experiments  examining  each  one  of 
the  above  points  are  presented  below: 

1.  The  strongest  correlation  is  shown  by  features 
computed  for  d  =  0.9375 ram,  which  is  the  interpixel 
distance.  To  see  therefore  whether  similar  results  could 
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Fig.  3.  Feature  |Ai,i|  in  grey  matter  for  d  =  0.9375 mm  versus 
the  result  of  the  MMSE  test  (•)  Alzheimer’s  Disease,  (+)  Control 
scans 

V.  Discussion 


MMSE 


be  obtained  by  performing  2D  analysis,  we  constructed 
the  2D  indicatrix  for  each  slice  separately  and  measured 
its  anisotropy: 

255  255 

M<M)  =  ;£^(fc-Z)2c(M;<M)  (8) 

k  1  1=1 

Figure  4  presents  an  example  of  a  2D  indicatrix. 


Alzheimer  Scan,  GDHL,  d=3mm 


x  axis 


In  this  section  we  are  trying  to  assess  the  validity  of 
our  approach  by  performing  various  tests.  In  particu- 


Fig.  4.  Example  of  a  2D  indicatrix 


TABLE  III 

The  best  features  identified  by  2D  analysis  on  Grey 

MATTER 


Best  Features 

Correlation 

with  MMSE 

l“4| 

d  = 

2 

-0.875 
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d  = 
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-0.859 
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-0.855 
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d  = 
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-0.828 
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d  = 
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0,Q 

d  = 

0.9375 

-0.821 

ag 

d  = 

0.9375 

-0.821 

ai 

d  = 

0.9375 

-0.817 

Features  F± ,  and  F2  can  be  computed  using  equations 
2  and  3,  and  F3  is  now  defined  by 


+  -my 

V  N 


(9) 


the  correlation  obtained  between  the  features  and  the 
MMSE  scores  was  not  greater  than  0.6. 

VI.  Conclusions 

The  method  of  graylevel  dependence  histograms 
(GLDH)  can  be  used  to  produce  many  features  strongly 
correlated  with  the  MMSE  scores  when  applied  to  the 
grey  matter  component  of  MRI-T1  scans. 

Our  experiments  showed  that  the  features  computed 
reflected  the  microtextural2  properties  of  the  brain 
rather  than  the  shape  of  each  region  of  interest,  and 
they  correlate  with  the  condition  of  the  patient  rather 
than  his/her  age. 

In  general  the  AD  brains  presented  greater 
anisotropy  in  their  grey  matter  texture  at  the  tested 
scales  than  control  brains.  2D  features  correlated  with 
the  MMSE  scores  as  well  as  3D  features,  pointing  that 
most  of  the  texture  information  is  already  available  in 
the  individual  layers. 

An  advantage  of  our  approach  is  that  only  one  scan 
is  required  to  compute  the  features.  This  is  a  differ¬ 
ent  approach  from  the  method  of  quantifying  volume 
change  between  registered  repeated  scans  as  a  diagnos¬ 
tic  marker  [4], 


Further  anisotropy  features  may  be  computed  from 
the  Fourier  descriptors  of  this  2D  shape: 

N/2 

R(<j>i)  =  a0  +  ( ak  cos (fc)  +  bk  sin(^))  (10) 

k=  1 

Table  III  presents  the  results  for  the  best  features 
computed  from  the  grey  matter  region  when  only  2D 
analysis  is  performed.  As  2D  analysis  is  performed  on 
each  slice  separately,  the  value  of  each  feature  was  aver¬ 
aged  over  all  slices  of  the  same  subject  before  attempt¬ 
ing  the  correlation  analysis. 

2.  The  correlation  between  the  texture  features 
and  the  age  of  the  subjects  (indepently  of  the  clini¬ 
cal  diagnosis)  is  smaller  than  0.6  (except  for  B2, 2  with 
correlation  of  0.674  at  d  =  0.9375 mm  and  0.694  at 
d  =  1.5 ram).  Furthermore  this  correlation  considering 
only  the  control  subjects  is  also  smaller  than  0.6  (ex¬ 
cept  for  A5i 2  with  correlation  0.626  at  d  =  0.9375 mm 
and  .85,1  with  correlation  of  0.613  at  d  =  3 mm).  It 
is  clear  that  for  these  particular  scans  the  process  of 
aging  was  not  strongly  reflected  in  the  texture  of  grey 
matter. 

3.  We  chose  the  distance  d  to  be  small  to  avoid  any 
interference  with  the  main  characteristics  of  the  brain 
(eg.  size,  shape  or  border  between  tissues).  In  order 
to  test  if  the  features  truly  measure  the  micro-textural 
characteristics  of  the  brain  rather  than  its  large  scale 
structure,  another  set  of  experiments  was  performed. 
The  values  of  the  voxels  of  the  grey  matter  were  sub¬ 
stituted  by  random  values,  and  both  texture  methods 
(3D-GLDH  and  2D-GLDH)  were  applied.  In  all  cases 
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